Why holes are not like electrons. III. How holes in the normal state turn into 

electrons in the superconducting state 
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In recent work, we discussed the difference between electrons and holes in energy band in solids 
from a many-particle point of view, originating in the electron-electron interaction[l[, and from a 
single particle point of view, originating in the electron-ion interaction^]. We proposed that super- 
conductivity in solids only occurs when the Fermi level is close to the top of a band (hole carriers), 
that it originates in 'undressing' of carriers from both the electron-electron and the electron-ion 
interaction, and that as a consequence holes in the normal state behave like electrons in the super- 
conducting stateQ]. However, the connection between both undressing effects was left unclear, as 
was left unclear how the transformation from hole behavior to electron behavior occurs. Here we 
clarify these questions by showing that the same electron-electron interaction physics that promotes 
pairing of hole carriers and undressing of carriers from the electron-electron interaction leads to 
undressing of carriers from the electron-ion interaction and transforms the behavior of carriers from 
hole-like to electron-like. Furthermore this phenomenon is connected with the expulsion of negative 
charge that we predict to occur in superconductors. These unexpected connections support the 
validity of our theoretical framework, the theory of hole superconductivity, to explain superconduc- 
tivity in solids. 

PACS numbers: 



I. INTRODUCTION 



When Felix Bloch formulated his epoch-making the- 
ory of electrons in metals he had an "uneasy feeling 
that the model of independent electrons might repre- 
sent a rather poor approximation and would turn out 
in some respects to be entirely inadequate" [f| (italic 
ours). Bloch realized full well the arbitrariness involved 
in his privileging the electron-ion interaction over the 
direct electron-electron Coulomb interaction, given that 
e 2 = lAAeVA is just as large an interaction between two 
electrons as between an electron and a monovalent ion 
at the same distance^. In the statement quoted above, 
with "in some respects" Bloch had undoubtedly the phe- 
nomenon of superconductivity in mind[6|, which was 
widely believed at the time to be caused by the electron- 
electron Coulomb interaction Q. However, Bloch's the- 
ory's myopic point of view became even more myopic 
with the conventional theory of superconductivity [8|], 
that adscribes the phenomenon to the interaction of elec- 
trons with ions displaced from their equilibrium position 
(electron-phonon interaction) Q , bringing the disregard 
for the role of the electron-electron Coulomb interaction 
to a new high. 

Instead, the theory of hole superconductivity [10] pro- 
poses that while electron-electron Coulomb interactions 
can indeed be neglected when a band is almost empty, 
they become increasingly dominant as the filling of a 
band increases. Unlike the prevalent point of view 
nowadays, that electron-electron interactions are most 
important near half- filling of a band 11| , we propose 
that electron-electron interaction physics dominates and 
qualitatively changes the normal metallic behavior when 



an electronic energy band is almost full, in particular giv- 
ing rise to superconductivity. 

We propose that the normal metallic and the supercon- 
ducting state are in a sense mirror images of each other, 
where the 'mirror' switches the sign of the electric charge, 
or the bottom and top of the band, or the electron-ion 
and the electron-electron interactions. The 'mirror' is by 
no means perfect because of the vastly different masses 
associated with the negative electron and the positive ion, 
hence there are significant differences between the normal 
metallic and the superconducting state, and the 'mirror 
switching' does not occur right at the half filled band but 
instead only when a band is almost full. Nevertheless, an 
essential commonality between both states emerges from 
considering the electron-ion and the electron-electron in- 
teraction on an equal footing, leading to the conclusion 
that normal metallic behavior ensues when the electron- 
ion interaction dominates, and superconducting behavior 
when the electron-electron interaction dominates. 

In the normal metal, there are few electrons in the band 
with a lot of room to move around, the positive ions are 
rigid and the electrons will act as nearly independent of 
each other and adjust their wavefunctions individually to 
optimize the electron-ion interaction, with the electron- 
electron interaction being non-optimized. In the super- 
conductor instead, the almost-full band is crowded with 
a lot of electrons that become 'rigid' because of their 
mutual strong interaction, with their wavefunction ad- 
justed to optimize the electron-electron interaction and 
with the electron-ion interaction being non-optimized. 
This means in particular that states that are electron- 
ion-interaction-energy costly will be occupied in the su- 
perconducting state. A consequence of this point of view 
is also that the theory predicts that superconductivity 
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is particularly favored when the electron-ion interaction 
stren gth Ze 2 is weak, i.e. when the ionic charge Z is 
small 

In the following, we analyze that aspect of the Coulomb 
interaction that we believe to be essential to understand 
superconductivity, at the level of a single atom, a di- 
atomic molecule and a solid. We show that the essen- 
tial physics manifests itself in both real and momentum 
space in a remarkably parallel fashion in going from the 
atom to the solid. It involves expansion of the electronic 
wavefunction to achieve lowering of electron-electron in- 
teraction and of quantum kinetic energy, at the expense 
of electron-ion interaction energy. 



II. THE ATOM 

Consider the wavefunction of an electron in the low- 
est energy state (Is) of a hydrogen-like atom of nuclear 
charge Z 
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(1) 



with r measured in units of the Bohr radius et . For two 
electrons in the atom in a singlet state, the spatial wave 
function is not 



(2) 



because the strong Coulomb interaction between two 
electrons in this state 



U = / d 6 rd 6 r'\ip z (n)\ 



In - r 2 \ 



-\^z(r 2 )\ 2 (3) 



makes this state energetically too costly. The true wave 
function in the two-electron atom contains both radial 
and angular correlations between the electrons, and is 
well described by the Hylleraas wave function How- 
ever, the main effect is captured by the Hartree approx- 
imation that allows for expansion of the single-particle 
orbital: 



with 



Z = Z-y 



(4a) 



(4b) 



The energy of the electrons in that state is (expressed in 
Rydbergs=13.6eV) 



with 



E(Z) = E km {Z) + Ee-i(Z) + E e _ e (Z) (5) 



E kin (Z) = 2Z 2 = E kln (Z) - AZy + 2y 2 (6a) 



E e _i(Z) = ~4ZZ = E e _i(Z) + AZy 



(6b) 
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FIG. 1: The two-electron atom. When the two electrons 
occupy the unexpanded orbital, they pay a high price in 
Coulomb energy, U. Instead, they prefer to reside in an ex- 
panded orbital, paying each the price e in single-particle en- 
ergy, yet achieving a lower total energy U 1 + 2e. 



E e _ e (Z) = - A Z = E^ e {Z)-^y (6c) 

with E e _ e (Z) = U. Therefore, expansion of the orbital 
costs electron-ion energy 

A£ e -, = 4Zy (7a) 

but gives a lowering of both kinetic energy 

AE kln = -AZy + 2y 2 (7b) 

and of Coulomb energy 
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(7c) 



The single-particle energy of each electron in the ex- 
panded orbital is the sum of its kinetic and electron-ion 
energy, and is larger than the single-particle energy in 
the unexpanded orbital by 



e= -(AE km + AE e ^) =y A 



(8) 



because the cost in electron-ion energy is larger than the 
kinetic-energy lowering. However, this is more than offset 
by the reduction in Coulomb repulsion in the expanded 
orbital: 



U' = E e _ e (Z) = U-~y 



so that the total energy 



U' + 2e = U + 2y 2 - -y 



(9) 



(10) 



is lower than the energy in the non-expanded doubly oc- 
cupied orbital provided y < 5/8. 

We represent the situation schematically in Fig. 1, 
with two "orbitals" per atom. The two electrons 
"choose" to occupy the orbital with higher single-particle 
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energy because the lower Coulomb repulsion more than 
offsets the cost in single-particle energy. The space- 
charge distribution becomes more extended in real space: 
negative charge is "expelled" outward when the second 
electron is added to the orbital, leading to the lowering 
of both electron-electron Coulomb energy (Eq. 7c)) and 
kinetic energy (Eq. (7b)). 

The optimal value of y that minimizes the energy 
Eq. (10) is y = 5/16, so the orbital Eq. (1) expands 
to (fz( r ) with Z = Z — y^, which costs electron-ion 
Coulomb energy AE e _i — jZ but gives a lowering of ki- 
netic energy AEkin = — + an d of electron-electron 
Coulomb energy AE e - e = — 1|. The single-particle en- 
ergy of each electron in the expanded orbital is larger 
than the single-particle energy in the unexpanded or- 
bital by e = ^7 which is more than offset by the re- 
duction in Coulomb repulsion in the expanded orbital 
U' = U(Z - 5/16) = U(Z) - §f, so that the total energy 



is U' + 2e = U- 



12S 



and the energy lowering achieved by 



expanding the orbital is 2b/128Ry = 2.7eV. 

The physics of the two-electron atom just described 
is a remarkable microcosm of the physics of electrons 
in a nearly filled electronic energy band: orbital expan- 
sion and promotion to higher single-particle energy levels 
driven by electron-electron repulsion and kinetic energy 
lowering will play a key role, as we discuss in the following 
sections. 



III. THE DIATOMIC MOLECULE 

Within a linear combination of atomic orbitals ap- 
proach, the bonding and antibonding orbitals of a di- 
atomic molecule are given by 



Vv( r ) = 



tpxjr) ± tp 2 {r) 
(2(1 ± S 12 ))V2 



(11) 



with the + (— ) sign corresponding to bonding (b) (anti- 
bonding (a)), assuming s-orbitals for definiteness. S±2 = 
(tpi,ip2) is the overlap matrix element. The bonding 
orbital has larger amplitude and correspondingly larger 
charge density in the region between the atoms, while 
the antibonding orbital changes sign and has vanishing 
charge density at a point between the atoms, as depicted 
in Fig. 2(a). As emphasized in II, this is a real physi- 
cal difference between bonding and antibonding orbitals 
that cannot be eliminated by a canonical transformation. 

Putting two electrons in the bonding orbital results in 
an electron-electron repulsion energy 



U bb = 



tFYdV|V>&(r)r— — ~ 
\r — r' 

U + V + 2J + AAt 



2(1 + S 12 )) 2 

and for two electrons in the antibonding orbital 
U + V + 2 J - 4Ai 



U a 



2(1 -S12)) 2 



(12a) 



(12b) 



diatomic molecule 




FIG. 2: (a) Schematic depiction of the wavefunction for bond- 
ing and antibonding orbitals in a diatomic molecule. The 
bonding orbital has higher charge density between the ions to 
get maximum benefit from the electron-ion interaction and a 
smooth wavefunction to give lower kinetic energy, (b) and (c) 
shows the bonding (lower) and antibonding (upper) single- 
particle energy levels, with different occupations of two elec- 
trons of opposite spin. 



where, in terms of the Coulomb integrals 



ripi(r>k(r) (13) 



(ij\kl) = / d A rd 6 r'<p*{r) V *{r') 



U = V = J = = and 

At = All these matrix elements are positive. In 

particular, the hybrid Coulomb matrix element 



At = (ii\ij) 



(14) 



lowers the Coulomb repulsion for two electrons in anti- 
bonding statesJUJ, as seen from Eq. (12b). For suffi- 
ciently large At, the 'inverted occupation' shown in Fig. 
2(c) with the two electrons in the antibonding orbital 
would have lower total energy (kinetic + electron-ion + 
electron-electron) than the usual one where both elec- 
trons are in the bonding orbital, Fig. 2(b). 

Note also that for one electron in the bonding state 
and one in the antibonding state, the direct Coulomb 
repulsion is 



U ba 



U + V-2J 



2(1 



2) 



(15) 



which indicates that the Coulomb matrix element J 
favors ferromagnetism, particularly near the half-filled 
band [3]. 

Let us now consider the diatomic molecule as a micro- 
cosm for an electronic energy band, as depicted in Fig. 
3. If the electron-ion energy dominates, the occupation 
will be as shown on the left side. In particular, for three 
electrons in the molecule (analogous to a 'nearly filled 
band') two electrons will go into the bonding state and 
one electron into the antibonding state. Instead, if the 
electron-electron interaction dominates, the occupation 
will be as depicted on the right side of Fig. 3, where 
two electrons occupy the antibonding state and one elec- 
tron the bonding state. Equivalently, the single hole in 
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given by Eq. (17): it is increasingly positive as the inter- 
atomic distance R decreases and the ionic charge Z de- 
creases. The 'correlated hopping' parameter At changes 
the hopping amplitude depending on the electronic oc- 
cupation of the sites involved in the hopping process: 
the hopping amplitude for an electron is to, to — At and 
to — 2At depending on whether there are 0, 1 or 2 other 
electrons at the sites involved in the hopping process. 

For one electron in the diatomic molecule, the eigen- 
states of the Hamiltonian Eq. (16) are the bonding and 
antibonding states c£ CT |0 >,c\ t7 \0 >, with |0 > the empty 
molecule and 



FIG. 3: The diatomic molecule as a microcosm for an energy 
band. When the electron-electron interaction dominates, the 
occupation of single-particle energy levels changes, as shown 
on the right side of the figure: the middle diagram represents 
a ferromagnet and the upper diagram a superconductor. 



the filled band resides in the bonding rather than in the 
antibonding state. 

In a tight binding description, the Hamiltonian for elec- 
trons in a diatomic molecule with one orbital per atom 
is given bv[l6| 

H = - }}tp - At(ni t - a + TO2,_o-)] + U/,niinji 

a i 

+ V n ^2a + J(^(c f 1(J C 2CT + h.C.)f (16) 
a a 

assuming higher order interactions (involving 6 or more 
fermion operators) can be neglected, to is the hopping 
amplitude for a single electron in the molecule. The cre- 
ation and annihilation operators in Eq. (16) refer to or- 
thogonal orbitals <j>\ (r) , (f>2 (r) obtained from linear com- 
binations of the atomic orbitals ld-\ (r), if2 (r)[16||. The 
Coulomb matrix elements are given by the expression 
Eq. (13) with the orthogonal orbitals. In particular, the 
hybrid matrix element At is given in terms of Eq. (13) 
with atomic orbitals, At a .o., by 



-1<T 



-2a 



At = At n 



Si 



■(U + V) 



(17) 



to lowest order in the nearest neighbor overlap 5*12. To 
the extent that the Mullikcn approximation for the over- 
lap charge distribution holds 17], the expression Eq. (17) 
exactly vanishes, and in practice it is found that devia- 
tions from the Mulliken approximation are very small [l6j. 
For that reason, the "correlated hopping" interaction 
term At is usually ignored in formulating tight binding 
Hamiltonians to describe interacting electrons in solids. 

However, in considering the interactions between elec- 
trons as the number of electrons in the diatomic molecule 
(or the band) increases it is essential to take into account 
the expansion of the atomic orbital that occurs for the 
doubly-occupied atom, discussed in Sect. II. When doing 
so, we showed in Refs. [TH, [ID, [l!| that the correlated 
hopping term At has a value very different from that 



^ba 



V2 



■■la 



-2a 



V2 



(18a) 



(18b) 



and energies ej, = —to, e a = +t . For three electrons in 
the molecule, the eigenstates are 



\b> c 



a >o 



t t 

C aa c b^bl 



l,|o> 



C at C al C ba 



|0> 



with energies 



h = -{t -2At) + U + 2V 



e a = +{t Q -2At) + U + 2V 



(19a) 



(19b) 



(20a) 



(20b) 



respectively. The ordering of these states depends on the 
magnitude of the single hole hopping amplitude 



t h =t - 2At 



(21) 



For th > 0, ib < e a and the state \b > CT , correspond- 
ing to two electrons in the bonding orbital and one elec- 
tron in the antibonding orbital has lower energy than the 
state \a > CT , so the lowest energy state corresponds to the 
three-electron state on the left side of Fig. 3. Instead, 
for th < the ordering is reversed and two electrons 
will occupy the antibonding orbital, as depicted on the 
right side of Fig. 3. The latter situation will occur if 
At > t /2. 

The parameter At represents a 'bond charge 
repulsion'[23| between electrons: it pushes electrons away 
from the bond (region between the ions) towards the sites 
(ions) by suppressing the occupation of the bonding state 
and increasing the occupation of the antibonding state. 
The magnitude of the hopping parameter to, which is de- 
termined by the strength of the electron-ion attraction, 
decreases as the ionic charge Z decreases. Instead, as 
already mentioned, the magnitude of At increases as 
the ionic charge decreases because of the increased or- 
bital expansion under double occupancy. So when the 
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electron-electron interaction dominates over the electron- 
ion attraction in the diatomic molecule, th changes sign 
and this leads to the inverted occupation shown on the 
right side of Fig. 3. 

Furthermore, as the electrons doubly-occupy the an- 
tibonding orbital, they have larger amplitude for their 
wavcfunction at the atomic site, rather than in the in- 
teratomic bond, as depicted by the upper picture in Fig. 
2(a). Larger on-site occupation implies a larger expan- 
sion of the atomic wavefunction, as depicted in the lower 
right panel of Fig. 1, and hence negative charge is ex- 
pelled outward, away from the interatomic region. This 
in turn also causes a relative lowering of kinetic energy, 
as discussed in Sect. II for the single atom. 

In summary, our analysis of the diatomic molecule 
shows strong parallels with the discussion in Sect. II for 
the single atom: when the electron-electron interaction 
dominates over the electron-ion interaction, the single- 
electron higher energy level (antibonding state) becomes 
more occupied and the single-electron lower energy level 
(bonding state) becomes less occupied. In other words, 
holes tend to increasingly occupy the lower-energy single- 
electron energy levels. At the same time, the negative 
charge spatial distribution increases outward, which im- 
plies a higher positive charge distribution in the region 
inside the structure. That is, negative charge is "ex- 
pelled" from the interior region of the molecule towards 
the exterior. The interaction matrix element responsible 
for this effect in the diatomic molecule (and in the solid as 
we will see)) is At, and the essential underlying physics is 
orbital expansion driven by electron-electron interaction 
(since, as mentioned, in the absence of orbital expansion 
the parameter At is found to be zero) and kinetic energy 
lowering (since the atomic orbital expansion gives kinetic 
energy lowering as seen in Sect. I). 



IV. FURTHER ANALYSIS OF At 

In this section we discuss in more detail the physical 
origin of the occupation inversion that we argue can oc- 
cur in the diatomic molecule, because the same principles 
will apply to the solid state. In Refs. [HI, EH, E3] we cal- 
culated from first principles the hopping amplitude for a 
single electron, to, and for a single hole, th, in a diatomic 
molecule, taking into account the orbital relaxation ef- 
fect. These hopping amplitudes were obtained from the 
difference in energy of an electron or a hole in the bond- 
ing and antibonding states. As can be seen for example 
in Fig. 7(c) of Ref. flU. t he hole hopping amplitude 
th (denoted by £2 in ref. [191]) goes to zero and in fact 
changes sign for sufficiently small interatomic distance R 
and small ionic charge Z. In this section we discuss in 
more detail how the change in sign of th comes about. 

In paper I of this series [l| we have emphasized the 
change in hopping amplitude with occupation that arises 
from the overlap matrix element of the expanded and un- 
expanded atomic orbital. For the Is orbital discussed in 



sect. II, the atomic overlap matrix element is given by 

{ZZfl 2 



S = {(pz,Vz) 



W)3 



(22) 



with Z = Z — 5/16. S is a decreasing function of the ionic 
charge Z, and as Z — > 5/16 it becomes arbitrarily small. 
The hopping amplitude for an electron when there are m 
other electrons in the two sites involved in the hopping 
process is t m = S m to, and has the same sign as to. If 
we define At as the difference in hopping amplitude for 
a single hole and when there is a second hole in one of 
the sites, it is given by 



At = t S(l - S) 



(23) 



This physics will promote pairing of hole carriers [2 lj. but 
will not change the sign of the hopping amplitude for 
holes, since for a single hole we have 



th = toS 



(24) 



of the same sign as to- 

However, there are additional contributions to the vari- 
ation of hopping amplitude with occupation, that arise 
from matrix elements of the Coulomb interaction involv- 
ing electrons on neighboring sites. In the absence of or- 
bital expansion these contributions essentially cancel, as 
discussed in Sect. Ill, however in combination with or- 
bital expansion they have a fundamental effect, and Eq. 
(24) ceases to be correct. We discussed in Ref.[12j a cal- 
culation of the hopping amplitude for a single hole in the 
diatomic molecule obtained by approximating the bond- 
ing and antibonding states of the three-electron diatomic 
molecule by the wavefunction 



4>b,a(ri,r2,r3) = 

<pi(ri)<pi(r2)<p2(r 3 ) ± (pi{ri)Cp 2 (r2)(p 2 {r 3 ) 
(2(1±5 12 ))V2 



(25) 



where (f>i{rj) and fi(rj) are the expanded and unex- 
panded orbitals for the j — th electron in atom i = 1, 2 
and the overlap is 



<5l2 = (ip 1 ,ip 1 ) 2 ((p 1 ,(p 2 )- 



(26) 



where (<pi,<pi) = S ((Eq. (22)). The energies of these 
states are 

(112|ff 3 |il2)±(112|g 3 |122) f97 s 

1 ± i5i2 

where the labels 1, 1 are shorthand for the orbitals <pi,<pi, 
etc. The three-electron Hamiltonian is given by 



H 3 = H 3sp + H 3ee 



(28a) 



with the single-particle and electron-electron parts given 
by 



H 3sp = hi + h 2 + h 3 



(28b) 
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(28c) 



where ft, is the single-particle Hamiltonian for the i-th 
electron and hij the interaction potential between the 
i-th and j-th electron, given by (in atomic units) 

2 2Z 2Z 

h = -V 2 

n r 2 



(29a) 



hu = — 



(29b) 



with ri,T2 the electronic coordinate relative to atom 1, 
2, and ri 2 the distance between electrons. The differ- 
ence between the energies Eq. (27) give the hole hopping 
amplitude th 



th 



(30) 



In the regime where the overlap S is small (small ionic 
charge Z) we can neglect the overlap Eq. (26), being 
proportional to S 2 . The difference in energies between 
the bonding and antibonding states is then determined 
by the second matrix element in Eq. (27), given by 

(112|H 3 |122) = (112|ff 3sp |122) + (112|ff 3ee |122) (31) 

The first term in Eq. (31) (single-particle term) is 



(112|ff 3sp |122) = (1, h, 2)S 2 + 2(1, h, l)SSi 



(32) 



and is negative, so it lowers the energy of the bonding 
state and raises the energy of the antibonding state. The 
first term in Eq. (32) corresponds essentially to Eq. (24), 
the second term would be absent if the orbitals at neigh- 
boring sites were constructed to be orthogonal to each 
other. 

The second term in Eq. (31) (electron-electron inter- 
action term) is 

(il2|ff 3ee |122) = 2(1, 1, h a , 1, 2)5 + (1, 2, h 12 , 1, 2)S H 

(33) 

and is positive, so it raises the energy of the bonding 
state and lowers the energy of the antibonding state. It 
is clear then that the state of inverted occupation (an- 
tibonding state) has lower electron-electron interaction 
energy than the ordinary state (bonding state) . The dif- 
ference in energy between the states is the hole hopping 
amplitude 



t h = -(l,h,2)S 2 -2{l,h,l)SSa 

- 2(l,T,/i 12 ,l,2)5-(l,2,/i 12 ,l,2)5i2 



(34) 



where the first two terms, involving the electron-ion in- 
teraction, are positive (including their sign), and the 
last two terms involving the electron-electron interaction 
are negative (including their sign) . When the electron- 
electron interaction dominates over the electron-ion in- 
teraction (for small ionic charge Z) the sign of th switches 




FIG. 4: Dependence of the sign of the single hole hopping 
amplitude in the diatomic molecule th on the ionic charge 
Z and the interatomic distance R for the Is orbital. In the 
region above the curve the single hole hopping amplitude has 
the same sign as the single electron hopping amplitude to 
(to > 0), corresponding to the ordering of energy levels shown 
on the left side of Fig. 3. In the region below the curve the 
sign of the single hole hopping amplitude is opposite to that 
of the single electron hopping amplitude and the occupation 
of the levels for three electrons in the diatomic molecule is 
inverted as shown on the right side of Fig. 3. 



from positive to negative, and the state of inverted occu- 
pation becomes the low-energy state. 

In Ref.[l2l] we only explored the parameter regime 
where th has the same sign as to, i.e. where there is no 
occupation inversion. For sufficiently small Z and small 
interatomic distance R, th is found to change sign, as 
shown in Fig. 4. (For the calculation in Fig. 4 the over- 
lap Eq. (26) was included hence there are other terms 
in th beyond those given in Eq. (34), see Ref.pjj]). The 
region th < corresponds to the electron-electron dom- 
inated regime shown on the right side of Fig. 3, where 
occupation inversion occurs. For further details on the 
form of the matrix elements entering in the calculation 
leading to Fig. 4 the reader is referred to the Appendix 
of Refill. 



V. THE SOLID 

We consider the Hamiltonian for electrons in a band 
in a tight binding representation given by [la] 



H = - (t - At(n i - a + nj- a ))(4 (T Cj a + h.c.) 

<ij>,cr 

+U 2^2 n it n il + ^ TliTij (35) 



<y> 



Performing a particle hole transformation 



h 



(36) 
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FIG. 5: The electronic states at the Fermi energy look very 
different when the Fermi level is near the bottom and near the 
top of a band in a normal metal. Near the bottom of the band 
the wavelength is large and the state is smooth, resembling 
the free-electron plane wave function, nearly unaffected by the 
electron-ion potential. Instead, near the top of the band the 
wavelength approaches a single lattice spacing and the state 
is strongly modified by the electron-ion interaction. 



the Hamiltonian becomes 
H = + (th + ^t(n i ,-a + n j ,- <T ))(hl T h j<T + h.c.) 

<ij>,a 

+U y~]nifriii + V ^2 n i n j ( 37 ) 



i <%]> 

where the number operators are now = h\ G hi a and 

t h = t a - 2At. (38) 

We have shown that this Hamiltonian leads to pairing 
of holes and to superconductivity in a BCS formalism in 
the regime of low hole concentration 21]. Furthermore 
we have discussed that the transition to superconductiv- 
ity driven by At is associated with "undressing" [22} : the 
effective mass of the carriers decreases and the quasipar- 
ticle weight increases, and a transfer of spectral weight 
from high to low frequencies takes place both in the single 
particle spectral function (detectable in photoemission 
experiments) and in the optical absorption spectrum [23|. 

In addition, in paper II of this series we have argued 
that superconductivity is also associated with 'undress- 
ing' of carriers from the electron- ion interaction, and that 
it involves a wavelength expansion: the wavelength of car- 
riers at the Fermi energy grows from a microscopic length 



interatomic distance) to a much larger 



wavelength. Here we explain how this occurs. 

Figure 5 shows schematically the form of the wavefunc- 
tion of the carriers at the Fermi energy when the band is 
almost empty and when it is almost full. We are assum- 
ing that the hopping amplitude to is positive, as occurs 
for s orbitals, which implies that the minimum in the 
band occurs for k — 0. When the Fermi level is close to 
the top of the band (hole carriers) the state is "bumpy" 
rather than smooth. As discussed in II, it is not a state 
well-suited to conduct electricity: it is strongly coupled 
to the lattice and is (nearly) Bragg-scattered by it, trans- 
ferring the momentum it acquires from an external field 
to the ionic lattice in the process, and acquiring momen- 
tum in opposite direction. It is highly "dressed" by the 
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FIG. 6: When the system goes superconducting, the anti- 
bonding states near the top of the band become fully occu- 
pied with electrons and the bonding states at the bottom of 
the band become empty. In other words, the holes "Bose- 
condense" to the bottom of the band. The promotion of 
electrons to the higher electron-ion energy states is driven 
principally by the electron-electron interaction. 



electron-ion interaction. In addition, as discussed in I, 
as the hole propagates it causes a large deformation in 
the atomic cloud where it lands (Fig. 1), and as a conse- 
quence it is also highly dressed by the electron-electron 
interaction. 

As the system goes superconducting the carriers at the 
Fermi energy undress from both the electron-electron [l| 
and the electron-ion interaction^ in order to conduct 
better. Their wavelength expands, and they no longer 
'see' the discrete electron-ion potential Figure 6 illus- 
trates how this occurs. 

Indeed, we propose that in the transition to supercon- 
ductivity, the electronic occupation of the band is shifted 
upward, to occupy all the single-particle energy levels up 
to the top of the band, leaving the lowest energy levels in 
the band empty. In other words, the holes that were oc- 
cupying the uppermost electronic levels "Bose-condense" 
and now occupy the lowest electron-ion energy levels. By 
so doing, the electronic states that are now at the bound- 
ary between occupied and empty states (dashed line on 
the right side of Fig. 6) become "smooth" states, of the 
'bonding type', with large amplitude in the region be- 
tween the ions(24|. 

How can such complete reorganization take place? 
Just like in the diatomic molecule, it will occur for 
low hole concentration if the hopping amplitude t(n) 
(n=number of electrons per site) changes sign: 



t(n) =t - nAt < 



(39) 



so that in particular the single hole hopping amplitude 
th = t - 2At < 0. Then, in the Hamiltonian Eq. (37) 
the lowest energy state for a hole occurs for k = rather 
than for k = ir. This, however, is not just a 'canonical 
transformation': the states near the bottom and the top 
of the band are qualitatively different. 

As function of temperature, the hopping amplitude 
t(n) can change sign if the number of electrons n in- 
creases as the temperature is lowered. Indeed we showed 
in Ref. [25[ that several anomalous properties of high 
T c cuprates are explained by the assumption that the 
number of holes in the Cu — O planes decreases as the 
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electron-electron dominated energy band = superconductor 




FIG. 7: In the superconducting state, the holes occupy 
smooth bonding states, that were previously occupied by elec- 
trons in the normal state. The holes 'ride' on top of the rigid 
negatively charged ions without disturbing their charge dis- 
tribution as they do in the normal state[2J]. The electron 
energy band depicted on the right side includes the electron- 
electron interaction energy nAt and therefore it is inverted. 



atom diatomic molecule solid 




FIG. 8: Single-electron energy levels (upper part of the figure) 
and real space charge distribution (lower part of the figure). 
As electrons move up the energy levels, negative charge moves 
outward in the atom, the diatomic molecule and the solid in 
the superconducting state. In the solid in the normal state, 
the charge distribution is spatially homogeneous. The charge 
redistribution is driven by electron-electron repulsion and ki- 
netic energy lowering, and is larger for smaller ionic charge Z 
in all three cases. 



temperature is lowered, and that such a charge transfer 
process between the planes and off-plane charge reser- 
voir atoms is aided by the Coulomb matrix element At. 
In other systems the change in n with temperature may 
be associated with redistribution of occupation between 
different bands without charge transfer between different 
regions in real space. 

In Figure 7 we represent the situation in an alternative 
way, where the 'energy band' now includes the electron- 
electron interaction energy nAt. Because the sign of the 
hopping amplitude has changed for the almost-full band, 
the lowest energy electronic states occur at k — it and 
are occupied by electrons, and the highest occur at k = 
and are unoccupied. Because the band is almost full, the 
carriers at the Fermi energy now have a smooth wave- 
function, as depicted on the left side of Fig. 7, just like 
the smooth wavefunction of carriers at the Fermi energy 
in the normal state of a metal when the band is almost 
empty. The holes ride 'on top' of the negative ions, with- 
out disturbing them and without noticing the discrete- 
ness of the potential. This corresponds to the 'mirror 
image' of the lower left panel of Fig. 5. 

This reorganization in energy level occupation has a 
counterpart in the real space charge distribution. Figure 
8 shows the situation schematically for the atom, the di- 
atomic molecule and the solid. As higher single-electron 
energy levels become occupied, electronic charge moves 
outward in the three cases. The theory of hole supercon- 
ductivity predicts that negative charge is expelled from 
the interior of the metal towards the surface as it un- 
dergoes the transition to the superconducting state[2?|. 
and that as a consequence an excess negative charge 
density exists within a London penetration depth of 
the surface of superconductors, with p_ = —H c /At{ or 
p_ = —H c i/4iT for type I and type II superconductors 
respectively [23] • 

The superconductor looks like a 'giant atom'[28| with 
expanded orbitals. 



VI. SUPERCONDUCTIVITY FROM HOLE 
UNDRESSING 

As discussed in II and Ref. [|[ , there is plenty of exper- 
imental evidence that dressed hole carriers in the normal 
state become undressed from the electron-ion interaction 
and behave like undressed electron carriers in the super- 
conducting state. How this occurs becomes clear from 
the considerations in this paper. 

Let us review the experimental evidence for undress- 
ing of carriers from the electron-ion interaction in the 
superconducting state: 

(1) Rotating superconductor. A superconducting body 
rotating with angular velocity cJ develops a uniform mag- 
netic field in the interior [2{|, given by 

2m e c _ 

B = —to (40) 

e 

with m e the free electron mass. This has been measured 
in conventional, heavy fermion, and high T c cuprate 
superconductors [30j. The magnetic field always points 
parallel to the angular velocity, as given by Eq. (40), 
never antiparallel, indicating that the superfluid carriers 
behave as negatively charged carners[3ll| . 

(2) Gyromagnetic effect: If a magnetic field B is sud- 
denly applied to a superconductor at rest, the entire body 
is found to start rotating with angular momentum 

-* m c — * 

Lbody = TT^VB (41) 

2ne 

with m e the free electron mass, and V the volume (32|. 
This angular momentum is always antiparallel to the 
applied magnetic field for the same reason as in (1) 
above [3lj. 

(3) Bernoulli potential: Where there is a spatial varia- 
tion of the superfluid velocity an electric field is expected 
to develop, the resulting potential is termed Bernoulli 
potential j33j. Experimental measurements (34j are con- 
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sistent with an electric field given by 



E = -V-m e vi 
e 2 8 



(42) 



where e and m e are the /ree electron charge and mass, 
and v s is the superfluid velocity. In particular the sign of 
the Bernoulli potential measured corresponds to negative 
charge carriers. 

(4) Hall effect: The Hall coefficient is found to be es- 
sentially always positive in the normal state, correspond- 
ing to hole carriers and to change its sign from pos- 
itive to negative at temperatures slightly below T c , in- 
dicating that carriers change from hole-like to electron- 
like[H. 

(5) Wavelength expansion: We have proposed an ex- 
planation for the Meissner effect in superconductors in- 
volving a wavelength expansion of the carriers at the 
Fermi energy 0, [37| from the microscopic length kf, 1 
to the mesoscopic length 2Xl (Xl ^London penetra- 
tion depth), corresponding to the wavevector change 
k F ^qo = 1/(2A L ). 

These experimental observations and theoretical con- 
siderations indicate that the carriers of electric current 
in the normal state, dressed hole carriers, morph into 
undressed electron-like carriers in the superconducting 
state. How does this happen? In II and Ref.0], we sug- 
gested that the antibonding electrons at the top of the 
Fermi distribution in the normal state condense to the 
bottom of the band in the transition to superconductiv- 
ity. However, this would conflict with the Pauli principle, 
because the states at the bottom of the band are occu- 
pied! 

Furthermore, it is important to point out that it would 
be incorrect to assume that all the electrons in the band 
become 'undressed' from the electron-ion interaction in 
the superconducting state. Because the superfluid weight 
n s that enters into the London penetration depth Xl 



1 



47m, 



(43) 



most definitely corresponds to the hole concentration 
rather than to the electron concentration in the band 
(which is much larger). It is essentially the same charge 
carrier density that carries the current in the normal 
state. 

So the holes in the normal state have to behave like 
electrons in the superconducting state. The solution of 
this puzzle is what was depicted in Fig. 6. In Fig. 9 we 
show the expected evolution with temperature of the oc- 
cupations in the band structure and in the Brillouin zone. 
As the metal is cooled into the superconducting state, the 
holes condense to the bottom of the band, which corre- 
sponds to long-wavelength states near the center of the 
Brilloin zone, which are smooth bonding states as de- 
picted on the left side of Fig. 7. So it may be said that 
the superfluid carriers are still holes, not electrons. How- 
ever, the sign of the effective mass for carriers at the 
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FIG. 9: Occupied states in the band (upper part of the fig- 
ure) and in the Brillouin zone (lower part of the figure) as 
the temperature is lowered. The shaded areas of the Brillouin 
zone are occupied by electrons. When the temperature is low- 
ered below T c a 'hole pocket' is pierced at k ~ that grows 
as the temperature is further lowered, while the hole pock- 
ets near the edges of the Brillouin zone become progressively 
smaller and eventually disappear at T = 0. At temperatures 
< T < T c the empty states near k ~ n correspond to the 
normal fluid hole carriers and the empty states near k ~ 
correspond to the superfluid hole carriers. 



bottom of the band is opposite to that of carriers at the 
top of the band. For that reason, these superfluid hole 
carriers behave like electrons, and in particular exhibit a 
negative Hall coefficient, as well as the other signatures 
of electron transport reviewed at the beginning of this 
section. 

Figure 9 also suggests a new interpretation of the 'two- 
fluid model' [38| of superconductivity. The normal quasi- 
particles are the holes near the top of the band, i.e. 
k ~ 7T, while the superfluid are the holes at the bottom 
of the band, with k ~ 0. The conventional BCS th eory 
is argued to be consistent with the two-fluid model [39j, 
however it does not provide a clear separation of both 
components as Fig. 9 does. 

Note also that the continuous process by which the 
holes move from the top to the bottom of the band as the 
superfluid condensate develops represents the momentum 
space counterpart of the orbit expansion discussed in con- 
nection with the Spin Meissner effect [371] . 



VII. TWO-ORBITAL MODEL 

In previous work we have introduced an electronic 
model with two orbitals per site to describe the essential 
physics of electron-hole asymmetry (electronic dynamic 
Hubbard model) [40] . The spacing between single-particle 
energy levels is e, and the interactions are such that two 
electrons will occupy predominantly the higher single- 
particle level because the Coulomb repulsion there (If) 
is much smaller than both in the lower level and between 
one electron in each level. This corresponds to the sit- 
uation depicted schematically on the right side of Fig. 
1. 
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FIG. 10: Propagation of a pair in the strong coupling limit. 
The fact that the holes undress completely and propagate as 
free electrons implies there is no relaxation of the charge cloud 
as depicted in (b), rather the propagation is as depicted in (a), 
with the holes occupying the lower atomic energy level. 



Consider a pair of holes propagating within this model 
in the limit where the size of the pair is a single lattice 
spacing. This corresponds to the strong coupling limit 
where the single hole hopping amplitude goes to zero as 
discussed in Ref . [ill] . In Fig. 10 we show two possi- 
ble ways the pair propagation can occur. As in Ref.floT] 
(Eq. (4)), we assume for simplicity that the bare hop- 
ping amplitude is the same between all nearest neigh- 
bor orbitals, t. If the propagation is as shown in Fig. 
10(b), the electrons in the singly occupied site relax to 
the lower orbital, and the hopping amplitude is reduced 
by the overlap matrix element S (Eq. 27(a) of Ref. [40]). 
Instead, if the propagation is as shown in Fig. 10(a), the 
hopping amplitude is not reduced. This corresponds to 
the case where the holes occupy the lowest orbital, as 
depicted on the right side of Fig. 6. 

The wavefunction for a pair of holes is a linear combi- 
nation of the states 



>i 



(44a) 



(44b) 



where i, j are nearest neighbor sites and z is the number 
of nearest neighbors to a site. Its energy is the lowest 
eigenvalue of the matrix 



#1 



or the matrix 



Ho 



-2<V^ 
-2tJz -U' 



_ -2tSyfz 
-2tS^/z -U'-2e 



for the case of Fig. 10(a) and 10(b) respectively. In 
the limit where S is very small, the condition on the 



parameters for the lowest eigenvalue of Hi to be smaller 
than the lowest eigenvalue of H 2 is 



eU' 



(45) 



Here, t is the bare electron hopping amplitude for an 
empty band. For a cubic lattice it is given by 



t = 



2m „a 2 



3.81eF 
a (A) 2 



(46) 



Therefore, the condition Eq. (45^ can be satisfied for rea- 
sonable parameters, e.g. a = 2A, t — 0.95eV, e — leV, 
U' = 5eV. In that case, the hole propagates without 
disturbing the background, as shown schematically in 
the left diagram of Fig. 7. The propagation shown in 
Fig. 10(b) corresponds to 'partial' undressing from the 
electron-electron interaction, while that shown in Fig. 
10(a) corresponds to full undressing 



VIII. DISCUSSION 

In this paper we have discussed the point of view that 
the competition between electron-ion interaction domi- 
nance versus electron-electron interaction dominance in 
solids leads to normal metallic or superconducting be- 
havior depending on whether the former or the latter 
wins. Electron-electron interaction will dominate when 
an electronic energy band is almost full, i.e. when the 
carriers in the normal state are hole-like. Correspond- 
ingly, superconductors in nature are found to have posi- 
tive Hall coefficient in the normal state[35] . Furthermore, 
the electron-electron interaction strength e 2 will domi- 
nate over the electron-ion interaction strength Ze 2 when 
the ionic charge Z is small. Correspondingly, high T c 
superconductivity is found in materials with highly neg- 
atively charged substructures (planes) containing nega- 
tive ions, like the (Cu — 2 ) = planes in the cuprates, the 
(Fe — As)~ planes in the arsenides and the B~ planes in 
MgB 2 . 

Here we restricted the discussion of the atom and the 
diatomic molecule to the simplest case of a Is orbital. 
However exactly the same physics should take place in 
bands originating from other atomic orbitals. For exam- 
ple we showed in Ref. [19( that the quantities of interest 
for 2p orbitals in a diatomic molecule behave very simi- 
larly to those in the Is orbital. With p orbitals, the roles 
of k ~ and k ~ tt switch, however the essential physics 
determined by whether a state is near the bottom or near 
the top of the band remains the same. 

When one talks about electronic energy bands one is 
implicitly privileging the electron-ion interaction over the 
electron-electron interaction. The states at the bottom 
of a band have low electron-ion energy, and those at the 
top of a band have high electron-ion energy. However, 
conversely, as we have argued in this paper, electrons re- 
siding in states near the bottom of the band have high 
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electron-electron repulsion energy, and those near the 
top of the band have low electron-electron repulsion en- 
ergy. It is only natural to assume that in certain cases 
a complete reorganization can occur and a new state 
will emerge that optimizes the electron-electron interac- 
tion instead of the electron-ion interaction. That, we 
propose, is the superconducting state: electrons redis- 
tribute their occupation in the band states to occupy the 
high electron-ion energy states, low electron-electron en- 
ergy states, and leave empty the low electron-ion energy 
states, high electron-electron energy states. The concept 
of "holes" provides a natural language to describe the 
process: holes condense to occupy the low electron-ion 
states, giving rise to hole superconductivity. 

The physics proposed in this paper also illustrates an 
even closer connection between the phenomena of su- 
perconductivity and metallic ferromagnetism than pre- 
viously suspected [42|. In our previous work we pro- 
posed that both phenomena originate in electronic "bond 
charge repulsion" [lU, [ill. H i , both lead to lowering of the 



carrier effective mass[23l |45| and "undressing" [22t |4 



and both are driven by off-diagonal matrix elements of 
the Coulomb interaction 47| , with ferromagnetism dom- 
inating near the half- filled band[l4[ and superconductiv- 
ity when the band is almost full. For metallic ferro- 
magnetism it was always clear that the phenomenon in- 
volves occupying some states that arc unfavorable for the 
electron-ion interaction and emptying some states that 
are favorable to the electron-ion interaction. Namely, the 
majority spin electrons occupy antibonding states that 
were empty in the non-ferromagnetic state, and some 
bonding states that were occupied by the minority spins 
in the normal state become empty in the ferromagnetic 
state[14{. With the interpretation of superconductivity 
proposed in this paper it is clear that something quite 
similar occurs in the superconducting state, with empty 
antibonding states becoming full and full bonding states 
becoming empty. 

In previous work we have emphasized the contribu- 
tion to At arising solely from modulation of the single 



particle hopping amplitude by the on-site orbital expan- 
sion overlap matrix elements (Frank-Condon factor) [l[. 
That phy sics is contained in the 'dynamic Hubbard 
model' [4^] with an on-site Coulomb repulsion modulated 
by a local boson degree of freedom, or in a purely elec- 
tronic model with two orbitals per site and only on-site 
interactions 40] . These models have much of the rele- 
vant physics: they give rise to pairing and superconduc- 
tivity driven by "undressing" from the electron-electron 



interaction 



n by "u 

Mm- 



However these models don't allow for 



a change in sign of the single particle hopping amplitude, 
thus will not lead to "undressing" from the electron- 
ion interaction. That physics requires inclusion of off- 
site Coulomb matrix elements in the Hamiltonian as dis- 
cussed in Sect. III. 

In summary, in this paper we have proposed a new 
physical picture to describe superconductors: that it is 
the holes that condense to the bottom of the electronic 
energy band when a system becomes superconducting. 
This naturally ties together several elements introduced 
earlier within the theory of hole superconductivity: it 
explains why carriers undress from the electron-electron 
and the electron-ion interactions and behave as com- 
pletely free electrons, yet their number is the number 
of holes rather than the number of electrons in the band; 
the physics is tied to the Coulomb matrix element At 
that h as p layed a key role since the beginnings of this 
theory [l5j ; finally, we have shown that this physics is con- 
nected to the negative charge expulsion from the interior 
of superconductors previously found based on different 
arguments^. 

In a finite cluster, this reorganization of energy level 
occupation should also occur, amongst the discrete en- 
ergy levels of the cluster. It has not escaped our notice 
that this may provide an explanation for the remarkable 
experimental observations of de Heer and coworkers [Ed] . 

Further discussion and development of this physics and 
its connection with other elements of the theory will be 
given in future work. 



[1] J.E. Hirsch, Phys.Rev. B65, 184502 (2002), heareafter [11] 

referred to as I. 
[2] J.E. Hirsch, Phys.Rev. B71, 104522 (2005) heareafter re- 
ferred to as II. [12] 
[3] J.E. Hirsch, Int. J. Mod. Phys. B 17, 3236 (2003). [13] 
[4] F. Bloch, Z. Phys. 52, 555 (1929). [14] 
[5] F. Bloch, Proc. R. Soc. London A 371, 24 (1980). [15] 
[6] F. Bloch, Handbuch der Radiologic, 2nd. ed., Vol. VI, p. [16] 
226 (1933). [17] 
[7] W. Meissner, Ergebnisse der exakten Naturwis- [18] 
senschaften 11, 219 (1932). [19] 
[8] J. Bardeen, L. N. Cooper, and J. R. Schrieffer,Phys. Rev. [20] 

108, 1175 (1957). 

[9] H. Froelich, Phys. Rev. 79, 845 (1950). [21] 
[10] See http://physics.ucsd.edu~jorge/hole.html for a list of 
references. 



"The Hubbard Model: Its Physics and Mathematical 
Physics", ed. by Dionys Baeriswyl et al. (Plenum, New 
York, 1995). 

J.E. Hirsch, Phys.Rev. B48, 3327 (1993). 

E. A. Hylleraas, Z. Phys. 54, 347 (1929) 

J.E. Hirsch, Phys.Rev. B40, 2354 (1989). 

J.E. Hirsch, Physica C 158, 326 (1989). 

J.E. Hirsch, Chem. Phys. Lett. 171, 161 (1990). 

R.S. Mulliken, J. Chim. Phys. 46, 497 (1949). 

J.E. Hirsch, Phys.Rev. B48, 3340 (1993). 

J.E. Hirsch, Phys.Rev. B48, 9815 (1993). 

S. Kivelson, W.-P. Su, J. R. Schrieffer, and A. J. Heeger, 

Phys. Rev. Lett. 58, 1899 (1987). 

J. E. Hirsch and F. Marsiglio, Phys.Rev. B 39, 11515 
(1989); F. Marsiglio and J.E. Hirsch, Phys. Rev. B 41, 
6435 (1990). 



12 



[22] J.E. Hirsch, Phys. Rev. B 62, 14498 (2000); Phys. Rev. 

B 62, 14487 (2000). 
[23] J.E. Hirsch and F. Marsiglio, Phys. Rev. B 62, 15131 

(2000); J.E. Hirsch, Physica C 199, 305 (1992). 
[24] J.E. Hirsch, Mat.Res.Soc.Symp.Proc. Vol. 156, 349 

(1989). 

[25] J.E. Hirsch and F. Marsiglio, Physica C 195, 355 (1992). 
[26] J.E. Hirsch, Phys.Rev.B 68, 184502 (2003). 
[27] J.E. Hirsch, Ann. Phys (Berlin) 17, 380 (2008). 
[28] J.E. Hirsch, Phys.Lett. A309, 457 (2003). 
[29] R. Becker, F. Sauter and C. Heller, Z. Phys. 85, 772 
(1933). 

[30] A.F. Hildebrand , Phys.Rev.Lett. 8, 190 (1964); N.F. 
Brickman, Phys.Rev. 184, 460 (1969); J. Tate, S.B. Felch 
and B. Cabrera, Phys.Rev. B42, 7885 (1990); A. A. Ver- 
heijen, J.M. van Ruitenbeek, R. de Bruyn Ouboter, and 
L.J. de Jongh, Physica B 165-166, 1181 (1990); Na- 
ture 345, 418 (1990); M.A. Sanzari, H.L. Cui and F. 
Karwacki, Appl. Phys. Lett. 68, 3802 (1996). 

[31] J.E. Hirsch, Phys.Rev.B 68, 012510 (2003). 

[32] I.K. Kikoin and S.W. Gubar, J. Phys. USSR 3, 333 
(1940); R.H. Pry. A.L. Lathrop and W.V. Houston, 
Phys.Rev. 86, 905 (1952); R. Doll, Zeits. f. Physik 153, 
207 (1958). 

[33] F. London, "Superfluids" , John Wiley & Sons, Inc., New 
York, 1950, Vol. 1. 

[34] C.J. Adkins and J.R. Waldram, Phys.Rev.Lett. 21, 76 
(1968); J. Bok and J. Klein, Phys.Rev.Lett. 20, 660 
(1968); T.D. Morris and J.B. Brown, Physica 55, 760 
(1971); Y.N. Chiang and O.G. Shevchenko, Low Temp. 
Phys. 22, 513 (1996). 

[35] I.Kikoin and B. Lasarev, Nature 129, 57 (1932); ZhETF 
3, 44 (1933); Physik.Zeits. d. Sowjetunion 3, 351 (1933); 



I.M. Chapnik, Sov,Phys. Doklady 6, 988 (1962). 
[36] H. Van Beelen, JP Van Braam Houckgeest, HM Thomas, 

C. Stolk, and R. De Bruyn Ouboter, Physica 36, 241 

(1967); C. H. Weijsenfeld, Phys.Lett. A 28. 362 (1968); 

N. Usui. T. Ogasawara, K. Yasukochi, and S. Tomoda, J. 

Phys. Soc. Japan 27, 574 (1969); K. Noto, S. Shinzawa 

and Y. Muto, Sol.St.Comm. 18, 1081 (1975); M. Galffy 

and E. Zirngiebl, Sol.St.Comm. 68, 929 (1988); S.J. Ha- 

gen, C. J. Lobb, R. L. Greene, M. G. Forrester and J. H. 

Kang, Phys.Rev.B 41, 11630 (1990); C.C. Almasan, S. 

H. Han, K. Yoshiara, M. Buchgeister, D. A. Gajewski, L. 

M. Paulius, J. Herrmann, M. B. Maple, A. P. Paulikas, 

Chun Gu, and B. W. Veal, Phys.Rev.B 51, 3981 (1995). 
[37] J.E. Hirsch, Europhys. Lett. 81, 67003 (2008); J. Phys. 

Cond. Matt. 20, 235233 (2008). 
[38] C.J. Gorter and H.G.B. Casimir, Phyz. Z. 35, 963 (1934). 
[39] J. Bardeen, Phys. Rev. Lett. 1, 399 (1958). 
[40] J.E. Hirsch, Phys.Rev. B 67, 035103 (2003). 
[41] J.E. Hirsch, Phyica C 161, 185 (1989). 
[42] J.E. Hirsch, Physica B 163, 291 (1990). 
[43] J.E. Hirsch, Phys.Rev.B 43, 705, Sect. V (1991). 
[44] D.K. Campbell. J.T. Gammel, and EY. Loh, Phys.Rev.B 

38, 12043, (1988). 
[45] J.E. Hirsch, Phys.Rev.B 59, 436 (1999). 
[46] J.E. Hirsch, Phys.Rev.B 62, 14131 (2000). 
[47] J.E. Hirsch, Physica C 341-348, 211 (2000). 
[48] J.E. Hirsch, Phys. Rev. Lett. 87, 206402 (2001). 
[49] J.E. Hirsch, Phys.Rev. B 65, 214510 (2002). 
[50] J.E. Hirsch, Phys.Rev. B 66, 064507 (2002). 
[51] S. Yin, S. Xu, A. Liang, J. Bowlan, R. Moro and W.A. de 

Heer, Jour, of Superconductivity and Novel Magnetism 

21, 265 (2008) and references therein. 



